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In this paper we compute the generating function of modular, fc-noncrossing diagrams. A 
fc-noncrossing diagram is called modular if it does not contains any isolated arcs and any arc has 
length at least four. Modular diagrams represent the deformation retracts of RNA pseudoknot 
structures (5J El [10] and their properties reflect basic features of these bio-molecules. The 
fNJ . particular case of modular noncrossing diagrams has been extensively studied [201 Ell 13 EJ- 

Let Qfc(n) denote the number of modular fc-noncrossing diagrams over n vertices. We derive 
exact enumeration results as well as the asymptotic formula Qfc(n) ~ Ckn~^- k ~ 1 ^ ~~^~l^ n f° r 
CN| ■ k = 3, . . . , 9 and derive a new proof of the formula Q 2 (n) ~ 1.4848 n" 3 / 2 1.8489" n ©. 
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1. Introduction 

X 

The main result of this paper is the computation of the generating function of modular 
fc-noncrossing diagrams, Q fe (z). A diagram is a labeled graph over the vertex set [n] = 
{1, ... ,n} with vertex degrees not greater than one. The standard representation of a 
diagram is derived by drawing its vertices in a horizontal line and its arcs in the 
upper half-plane. A fc-crossing is a set of k distinct arcs ji), (12,32), ■ ■ ■ , (ik,jk) with 
the property 

fi < %2 < ■ ■ ■ < ik < ji < 32 < ■ ■ ■ < 3k- 
Similarly a /c-nesting is a set k distinct such arcs such that 

z'i < 12 < . . . < ik < 3k < ■ ■ ■ h < 3i- 

Let A, B be two sets of arcs, then A is nested in B if any element of A is nested in any 
element of B. A diagram without any /c-crossings is called a /c-noncrossing diagram. The 
length of an arc, is s = j — i, and we refer to such arc as s-arc. Furthermore, 
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• a stack of length a, S?j, is a maximal sequence of "parallel" arcs. A diagram, in 
which any arc is contained in a s-stack, where s > 2 is called a canonical diagram, 

(i + 1, j - 1), . . . , (z + (a - 1), j - (a - 1))). 

S"j is also referred to as a a-stack. 

• a stem of size s is a sequence 

( CO"l o°"2 q<T a \ 

\ h Jl' »2j2' ' ' ' ' °t a ,jj 

where S'f" 1 ,- is nested in S^ m ~ 1 4 such that any arc nested in Sf m ~ X : is either 
contained or nested in S7 m „- , for 2 < m < s. 

• a hairpin-loop (loop) consists of an arc and a sequence of consecutive isolated 
vertices {z + 1, i + 2, . . . ,j — 1}, see Fig. CD 




3' 5 10 15 20 25 5' 

Fig. 1. Features of a modular 3-noncrossing diagram represented as planar 
graph (top) and in standard representation (bottom). We display a stack of 



length two (green), a stem of size two (red) and a hairpin- loop (purple). 

RNA secondary structures [201 [2H El [22] are in the language of diagrams exactly mod- 
ular, 2-noncrossing diagrams. In [H [91 [TOj [12] , various classes of /c-noncrossing diagrams 
have been enumerated. However the approach employed in these papers does not work for 
modular /c-noncrossing diagrams. In contrast to the situation for RNA secondary struc- 
tures, the combination of minimum arc length and canonicity poses serious difficulties. 
The main idea is to build modular /c-noncrossing diagrams via inflating their colored, 
Vfc-shapes, see Fig. [21 These shapes will be discussed in detail in Section [H The inflation 
gives rise to "stem-modules" over shape-arcs and is the key for the symbolic derivation 
of Qfc(z). One additional observation maybe worth to be pointed out: the computation 
of the generating function of colored shapes in Section HI hinges on the intuition that the 
crossings of short arcs are relatively simple and give rise to manageable recursions. The 
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coloring of these shapes then allows to identify the arc-configurations that require special 
attention during the inflation process. Our results are of importance in the context of 
RNA pseudoknot structures [15] and evolutionary optimization [?]. Furthermore they 
allow for conceptual proofs of the results in [U |9l [101 [T2] . 





Fig. 2. Modular /c-noncrossing diagrams: the inflation method. A modular 3- 
noncrossing diagram (top) is derived by inflating its V3-shape (bottom) in two 
steps. First we individually inflate each shape- arc into a more complex configu- 
ration and second insert isolated vertices (purple). 



The paper is organized as follows. In Section [2] we recall some basic facts on singularity 
analysis, the generating function of fc-noncrossing matchings, V^-shapes and symbolic 
enumeration. Since the results on /c-noncrossing matchings are difficult to recover from 
the literature or simply new but immaterial for our purposes, we present their proofs in 
the Supplemental Materials (SM). In Section [3] we analyze modular, noncrossing diagrams 
and in Section H] we compute the generating function of colored shapes. Finally we prove 
the main theorem in Section [5j 
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2. Some basic facts 



2.1. Singularity analysis. Oftentimes, we are given a generating function without hav- 
ing an explicit formula of its coefficients. Singularity analysis is a framework that allows 
to analyze the asymptotics of these coefficients. The key to the asymptotics of the co- 
efficients are the singularities, which raises the question on how to locate them. In the 
particular case of power series f(z) = J2 n>0 a n z n with nonnegative coefficients and a 
radius of convergence R > 0, a theorem of Pringsheim [2J [T7], guarantess a positive real 
dominant singularity at z — R. As we are dealing here with combinatorial generating 
functions we always have this dominant singularity. We shall prove that for all our gener- 
ating functions it is the unique dominant singularity. The class of theorems that deal with 
the deduction of information about coefficients from the generating function are called 
transfer-theorems [2]. 



Theorem 1. [2] Let [z n ]f(z) denote the n-th coefficient of the power series f(z) at z = 0. 
(a) Suppose f\z) = (1 — z)~ a , a G C \ Z< , then 

' a_1 r a{a-l) a{a - l)(a - 2) (3a - 1) 
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rial 



1 + 



2n 
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(2.1) 



a 2 (a - l) 2 (a - 2) (a -3) Q fl 



n' 



(b) Suppose f(z) = (1 — z) r \og(j^), r e Z> , then we have 
(2-2) [z n }f(z) ~ (-1)' 



n(n — 1) ... (n — r) 



We use the notation 

(2.3) (f(z) = 6 (g(z)) as z -> p) ^ (f(z)/g(z) ^cas^p), 

where c is some constant. We say a function /(z) is A p analytic at its dominant singularity 
z — p, if it analytic in some domain A p (<f), r) = {z \ \z\ < r,z ^ r, |Arg(^ — p)| > 0}, for 
some <p, r, where r > |p| and < < ^. Since the Taylor coefficients have the property 

(2.4) V 7 GC\0; [z n ]f(z)= 1 n [z n ]f(-), 

7 

We can, without loss of generality, reduce our analysis to the case where z — 1 is the unique 
dominant singularity. The next theorem transfers the asymptotic expansion of a function 
near its unique dominant singularity to the asymptotic of the function's coefficients. 
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Theorem 2. j2] Let f(z) be a Ai analytic function at its unique dominant singularity 
z = 1. Let 



1 - z ( 

That is we have in the intersection of a neighborhood of 1 

(2.5) f(z) = Q(g(z)) forz^l. 
Then we have 

(2.6) [z n ]f(z) = ([z n ]g(z)) . 

2.2. &;-noncrossing matchings. A /c-noncrossing matching is a /c-noncrossing diagram 
without isolated points. Let f k (2n) denote the number of fc-noncrossing matchings. The 
exponential generating function of /c-noncrossing matchings satisfies the following identity 

- = det[/ i _ i (2z)-/ i+J -(2^)] li] 



(2.7) Y^fk^n)-^— = det[I^(2z) - I i+j (2zr' 



n>0 

where I r (2z) = Ylj>o j\(j+r)\ * s ^ ne hyperbolic Bessel function of the first kind of order 
r. Eq. (12 .7\\ combined with the fact that recursions for the coefficients of the exponential 
generating function translate into recursions for the coefficients of the ordinary generating 
function, allows us to prove: 



Lemma 1. The generating function of k-noncrossing matchings over2n vertices, F k (z) = 
^ n>0 /fc(2n) z n is D-finite, [IB] , i.e. there exists some e G N such that 

(2.8) q^ k {z)^F k {z) + q lth ( z )±^Y k {z) +■■■ + q^ k (z)¥ k {z) = 0, 

where qj, k {z) are polynomials. 



We sketch the proof of Lemma Q] in the SM. 

Lemma [1] is of importance for two reasons: first any singularity of F k (z) is contained 
in the set of roots of qo,k(z) [16J, which we denote by R k . Second, the specific form of the 
ODE in eq. (12 .8p is the key to derive the singular expansion of F k (z), see Proposition [T] 
below. 

We proceed by computing for 2 < k < 9, the polynomials qo tk (z) and their roots, see 
Table [T] and observe that [TT] 

(2.9) f k (2n) ~ c k n -((k-±) 2 +(k~i)/2) ( 2 (k _ 1)^ z k > 0, k > 2. 
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k qo,k(z) R, 



2 
3 
4 
5 
6 
7 
8 
9 



lUz 2 -40z + l)z 3 {J, ^} 

1024^-80^ + 1)^ {i^} 
14400z 3 - 41442 2 + 1402 - l)z 5 M> Ton) 

1474562 3 - 12544z 2 + 2242 - l)z 6 {i i -^j 

28224002 4 - 8266242 s + 315842 2 - 3362 + 1)/ {* ^ ^ 

377487362 4 - 33587202 3 + 698882 2 - 4802 + 1)2 8 {1, ^, j^} 



TABLE 1. We present the polynomials %k{ z ) an d their nonzero roots obtained 
by the MAPLE package GFUN. 



Equation ( 12. 9 p guarantees that F k (z) has the unique dominant singularity pi, where 
Pk = l/(2k — 2). According to Lemma [IJ F k (z) is D-finite, whence we have analytic 
continuation in any simply connected domain containing zero avoiding its singularities 
[T9] . As a result Fk(z) is A p 2 analytic as required by Theorem [2j Lemma [U and eq. ( 12. 9 h 
put us in position to present the singular expansion of F k (z): 

Proposition 1. For 2 < k < 9, the singular expansion ofF k (z) for z — >• p\ is given by 

F (z) = [Pk{z - Pi) + C k (z - log(, - p2) (1 + o(1)) 

\ A(* - pD + <4(* - pl)^-^-^ (i + (i)) , 

depending on k being odd or even. Furthermore, the terms P k {z) are polynomials of degree 
not larger than (k — l) 2 + (k — l)/2 — 1, c' k is some constant, and pk = I /(2k — 2). 

We give the proof of Proposition [T] is the SM. 

In our derivations the following instance of the supercritical paradigm [2] is of central 
importance: we are given a D-finite function, f(z) and an algebraic function g(u) satis- 
fying g(0) = 0. Furthermore we suppose that f(g{u)) has a unique real valued dominant 
singularity 7 and g is regular in a disc with radius slightly larger than 7. Then the su- 
percritical paradigm stipulates that the subexponential factors of f(g(u)) at u — 0, given 
that g{u) satisfies certain conditions, coincide with those of f(z). 

Theorem [TJ Theorem [2] and Proposition [1] allow under certain conditions to obtain the 
asymptotics of the coefficients of supercritical compositions of the "outer" function F k (z) 
and "inner" function ip(z). 

Proposition 2. Let ip(z) be an algebraic, analytic function in a domain D = {z\\z\ < 
r} such that ip(0) = 0. Suppose 7 is the unique dominant singularity of ' F k (ijj(z)) and 
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minimum positive real solution 0/^(7) = p\, I7I < r, where ^'(j) 7^ 0. Then Fk{ijj(z)) 
has a singular expansion and 

(2.10) [z n ]F k (^(z)) ~ ^ n -((fc-i) a +(*-i)/2) , 

where A is some constant. 

2.3. Shapes. 

Definition 1. A VVshape is a fc-noncrossing matching having stacks of length exactly 
one. 

In the following we refer to V^-shape simply as shapes. That is, given a modular, k- 
noncrossing diagram, 5, its shape is obtained by first replacing each stem by an arc and 
then removing all isolated vertices, see Fig. [3] 




1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 1 2 3 4 5 6 



"A 




1 2 3 4 5 6 7 



Fig. 3. From diagrams to shapes: A modular, 3-noncrossingd diagram (top-left) 
is mapped in two steps into its vVshape (top-right). A stem (blue) is replaced 
by an single shape-arc (blue). 



Let Jfc(s,m) (ik(s,m)) denote the set (number) of the V^-shapes with s arcs and m 
1-arcs having the bivariate generating function 



,m)z s u m . 



(2.11) I fc (z,u) = ^^i fc (s, 

s>0 m=0 

The bivariate generating function of ik{s,m) and the generating function of F k (z) are 
related as follows: 

Lemma 2. [H] Let k be natural number where k > 2, then the generating function Ik(z, u) 
satisfy 

, 01 m y / \ 1 + z ( z(l + z) 

( 2 - 12 ) h{z,u) = — — F k 



l + 2z-zu \{l + 2z-zu) 
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2.4. Symbolic enumeration. In the following we will compute the generating functions 
via the symbolic enumeration method j2]. For this purpose we need the notion of a 
combinatorial class. A combinatorial class (C, Wg) is a set together with a size-function, 
wq\ C — > Z + such that C n = (n) is finite for any n G Z + . We write w instead of wq 
and set C n = |C n |. Two special combinatorial classes are £ and Z which contain only one 
element of size and 1, respectively. The generating function of a combinatorial class 6 
is given by 

(2.13) C(z) = J2 zW(c) = J2 C - Z ^ 

cee n>o 

where C n C 6. In particular, the generating functions of the classes £ and Z are E(z) = 1 
and Z(z) = z. Suppose C, D are combinatorial classes. Then C is isomorphic to D, 
6 = D, if and only if \/ n > 0,|C n | = \T> n \. In the following we shall identify isomorphic 
combinatorial classes and write C = T> if C = D. We set 



(2.14) w e+ v(a) 



C + D := 6 U D, if C n D = and for a G C + D, 

we(a) if a G C 
wd(«) if a G D. 
• 6 x D := {a = (c, d) \ c G C, d G T>} and for a G C x D, 
(2.15) ^exn((c, d)) = w; e (c) + 

and furthermore G m := H™ =1 Q and Seq(C) := £ + C + C 2 + • • • . Plainly, Seq(C) is 
a combinatorial class if and only if there is no element in C of size 0. We immediately 
observe 

Proposition 3. Suppose A, Q and D are combinatorial classes with generating functions 
A(z), C(z) and D(z). Then 

(a) A = Q + 'D^ A{z) = C(z) + B(z) 

(b) A = C x D =^ A(z) = C(z) ■ T>(z) 

(c) A — Seq(C) A(z) = iz^iy- 

3. Modular, noncrossing diagrams 

Let us begin by studying first the case k = 2 [6], where the asymptotic formula 

Q 2 (n) ~ 1.4848 • n" 3/2 • 1.8489 n , 

has been derived. In the following we extend the result in [6] by computing the generat- 
ing function explicitely. The above asymptotic formula follows then easily by means of 
singularity analysis. 



Proposition 4. The generating function of modular, noncrossing diagrams is given by 

, s n , . l-z 2 + z 4 ( z 4 -z 6 + z s 

(3-1) Q*{z) = , : ~ 2 , _ 3 , ^ ; _ B • F 



1 - z - z 2 + z 3 + 2z* + 2 6 V(l - z - 22 + z 3 + 2 ^ 4 + z 6 ) 5 
and iae coefficients Q2(n) satisfy 

Q 2 (n) ~ c 2 n- 3 / 2 72 -", 
where 72 i/ie minimal, positive real solution of $(z) = 1/4, and 

(3 ' 2) ^ = (l- 2 -^ + 2 3 + 2 ^ + 2 6)2- 

i/ere we /iave 72 ~ 1.8489 and c 2 ~ 1.4848. 



Proof. Let Q2 denote the set of modular noncrossing diagrams, J2 the set of all V2-shapes 
and ?2( m ) those having exactly m 1-arcs. Then we have the surjective map 

ip: Q 2 -> %. 

The map tp is obviously surjective, inducing the partition Q 2 = U 7 <^ -1 (7), where <£> _1 (7) 
is the preimage set of shape 7 under the map ip. Accordingly, we arrive at 

(3.3) Q2(z) = J2 E QtW- 

m>0 7 e J 2 (m) 

We proceed by computing the generating function Q 7 (z). We shall construct Q 7 (z) from 
certain combinatorial classes as "building blocks". The latter are: M (stems), % (stacks), 
INT (induced stacks), L (isolated vertices), "R (arcs) and Z (vertices), where Z(z) = z and 
R(z) = z 2 . We inflate 7 G $2(171) having s arcs, where s > max{l,m}, to a modular 
noncrossing diagram in two steps. 

Step I: we inflate any shape-arc to a stack of length at least 2 and subsequently add 
additional stacks. The latter are called induced stacks and have to be separated by means 
of inserting isolated vertices, see Fig. HJ Note that during this first inflation step no 
intervals of isolated vertices, other than those necessary for separating the nested stacks 
are inserted. We generate 



sequences of isolated vertices £ = Seq(Z), where 

1 



L(z) 



1 - z 
stacks, i.e. 

% = ti 2 x Seq (ft) 
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n ^ (1) 



1234 12345678 ^ 



1 2 3 4 5 6 7 8 9 10 11 12 13 




(2) 

1 2 3 4 5 6 7 8 9 10 11 12 13 



— (3) 

1 2 3 4 5 6 7 8 9 10 11 12 13 14 

Fig. 4. A shape (left) is inflated to a noncrossing, canonical diagram. 



with the generating function 



1 



induced stacks, i.e. stacks together with at least one nonempty interval of isolated 
vertices on either or both its sides. 

N = X x (ZxL + ZxL + (Zx L) 2 ) 

with generating function 

4 / / \ 2 N 

„w n z z l z 

N z) = -. ^ 2- 



1 — z 2 \ 1 — z ' \1 — z 

stems, that is pairs consisting of a stack % and an arbitrarily long sequence of 
induced stacks 

M = % x Seq (N) 

with generating function 

K(z) 



M(z) 



l-z 2 



Step II: we insert additional isolated vertices at the remaining (2s+l) positions. For each 
1-arc at least three such isolated vertices are necessarily inserted, see Fig. [51 Combining 
Step I and Step II we arrive at 

(3.4) Q 7 = (M) s x x (Z 3 x L) m , 
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1 2 3 4 5 6 7 8 910111213 1 2 3 4 5 6 7 8 910111213141516171819 20 

Fig. 5. Step II: the noncrossing canonical diagram (left) obtained in (1) in Fig.H] 
is inflated to a modular noncrossing diagram (right) by adding isolated vertices 
(red). 



where Q 7 is the combinatorial class of modular noncrossing diagrams having shape 7. 
Therefore, 

s 



Q 7 (z) 



1=JI ( 2 T3I + (1=;) ) 



1 \ 2s+l — m. / ■? 

1 \ / z d 



1 - z7 \l-z 



1 (l-z 2 )(l-z) 2 -(2z-z 2 )z 4 ' 



Z 4 



Since for any 7,7! £ 3 2 (s,m) we have Q 7 (z) = Q 71 (z), we derive 

m>0 7 6J 2 (m) s>0 m=0 

We set 

V ^ {l-z 2 )(l-z) 2 -{2z-z 2 )z i 
and note that Lemma [2] guarantees 

YY i 2 (s,m)x s y m = 1 + X Yf 2 (2s,0)( , ^ + %) , 9 Y 
Therefore, setting x = r)(z) and y = z 3 we arrive at 

_ , N 1 - 2 2 + ^ 4 / 2 4 - 2 6 + £ 8 



1 - z - z 2 + z 3 + 2z A + z 6 V (! - 2 - ^ 2 + z3 + 2zA + ^ 6 ) 2 
By Lemma [U Q 2 (z) is D-fmite. Pringsheim's Theorem [17] guarantees that Q, 2 {z) has a 
dominant real positive singularity 72. We verify that 72 which is the unique solution of 
minimum modulus of the equation $(z) = p 2 , where p\ is the unique dominant singularity 
of F 2 (z) and p 2 = 1/2. Furthermore we observe that 72 is the unique dominant singularity 
of Q 2 (z), see the SM. It is straightforward to verify that $'(72) 7^ 0. According to 
Proposition |2} we therefore have 

Q 2 (n) ~ c 2 n- 3 / 2 7 2 " n , 

and the proof of Proposition H] is complete. □ 
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4. Colored shapes 

In the following we shall assume that k > 2, unless stated otherwise. The key to com- 
pute the generating function of modular fc-noncrossing diagrams are certain refinements 
of their VVshapes. These refined shapes are called colored shapes and obtained by dis- 
tinguishing a variety of crossings of 2-arcs, i.e. arcs of the form (i, i + 2). Each such class 
requires its specific inflation-procedure in Theorem [3j 

Let us next have a closer look at these combinatorial classes (colors): 

• Ci the class of of 1-arcs, 

• C 2 the class of arc-pairs consisting of mutually crossing 2-arcs, 

• C3 the class of arc-pairs (a, (3) where a is the unique 2-arc crossing f3 and f3 has 
length at least three. 

• C4 the class of arc-triples (ai,(3,a 2 ), where a.\ and ot 2 are 2-arcs that cross j3. 

In Fig. [6] we illustrate how these classes are induced by modular fc-noncrossing diagrams. 




Fig. 6. Colored vVshapes: a modular 3-noncrossing diagram (top) and its col- 
ored VVshape (bottom). In the resulting vVshape we color the four classes as 
follows: Ci(green), C2(black), C3(blue) and C^red). 



Let us refine V fc -shapes in two stages. For this purpose let Jfc(s, Mi, w 2 ) and ik(s,Ui,u 2 ) 
denote the set and cardinality of V^-shapes having s arcs, U\ 1-arcs and u 2 pairs of 
mutually crossing 2-arcs. Our first objective consists in computing the generating function 

s L^J 

W k (x,y,w) = Yl Yl *k(s,u u u 2 )x s y Ul w U2 . 

s>0 ui=0 U2=0 

That is, we first take the classes Ci and C 2 into account. 
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Lemma 3. Fork > 2, the coefficients ik(s,Ui,U2) satisfy 

(4.1) i k (s,u 1 ,u 2 ) = foru 1 + 2u 2 >s 

(4.2) ^ ik{s,u 1 ,u 2 ) = i k (s,Ui), 

112=0 

where i k (s,ui) denotes the number of\l k -shapes having s arcs, u\ 1-arcs. Furthermore we 
have the recursion: 

(u 2 + l)ik(s + l,u 1 ,u 2 + 1) = (ui + l)i fc (s,ui + 1, u 2 ) 

(4.3) +(ui + l)i fc (a - 1, mi + 1, m 2 ). 

The proof of Lemma [3] is given in SM. We next proceed by computing ~W k (x,y,w). 
Proposition 5. For k > 2, we have 

(4.4) W fc (x, y, w) = (1 + a> F fc (x(l + x)t> 2 ) , 
where v — ((1 — w)x 3 + (1 — w)x 2 + (2 — y)x + 1) 1 . 

Proof. According to Lemma [2} we have 

h(z,u) = — — F fc ' 



l + 2z-zu~" \{l + 2z-zu) 2 J ' 

This generating function is connected to W/ C (x, y,z) via eq. (14. 2 p of Lemma |3] as follows: 
setting w = 1, we have W k (x,y, 1) = Ifc(x,y). The recursion of eq. (14. 3 p gives rise to the 
partial differential equation 

^ 45 ^ dW k (x,y,w) _ ^ dW k (x,y,w) | ^ 2 <9W fc (x, w) 



<9u> <% <% 

We next show 

• the function 

(1 + x) 



W* k (x,y,w) 

(4.6) 



(1 — w)x 3 + (1 — w)x 2 + (2 — y)x + 1 
(1 + x)x 



((1 - w)x 3 + (1 - w)x 2 + (2 - y)x + If 
is a solution of eq. (14. 5p . 

its coefficients, i* k (s,Ui,u 2 ) = [x s y Vl w V2 ]Wl(x, y, w), satisfy 

i k (s, ui, u 2 ) — for u\ + 2u 2 > s, 
Wl(x,y,l) = l k (x,y). 
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Firstly, 

(4.7) ™i(x,y,w) = uFk{u) + 2uK{u) 

(4.8) dW* k (x,y,w) = + ^ + + , 

ow 



where 



x(l + x 

u 



((1 - w)x 3 + (1 - w)x 2 + (2 - y)x + l) 2 



and F' k (u) = J2n>o n fk(2n)(u) n . Consequently, we derive 
^ 4g ^ dW* k {x,y,w) _ x dWj{x,y,w) | ^ 2 dW* k {x,y,w) 



dw dy dy 

Secondly we prove i* k (s,ui,u 2 ) = for n x + 2u 2 > s. To this end we observe that 
W£(x, y,w) is a power series, since it is analytic in (0,0,0). It now suffices to note that 
the indeterminants y and w only appear in form of products xy and x 2 w or x 3 w. Thirdly, 
the equality W£(x, y, 1) = l k (x, y) is obvious. 
Claim. 

(4.10) W* k (x,y,w)=W k (x,y,w). 

By construction the coefficients i£(s, u\, u 2 ) satisfy eq. (14.31) and we just proved il(s, ui, u 2 ) - 
for u\ + 2u 2 > s. In view of Wj£(a;, y, 1) = Ifc(x, y) we have 

Vs,«i; ^ i£(s,Ui,u 2 ) = «fc(s,«i). 

M2=0 

Using these three properties it follows via induction over s 

Vs,u 1 ,u 2 > 0; il(s,Ui,u 2 ) = ik(s, u u u 2 ), 
whence the Claim and the proposition follows. □ 



In addition to Ci and C 2 , we consider next the classes C 3 and C 4 . For this purpose 
we have to identify two new recursions, see Lemma III Setting u = (ui, . . . , u 4 ) we denote 
by "J k (s,u) and i k (s,u), the set and number of colored V^-shapes over s arcs, containing 
Ui elements of class Cj, where 1 < % < 4. The key result is 

Lemma 4. For k > 2, the coefficients i k (s,u) satisfy 

(4.11) i k (s,Ui,u 2 ,u 3 ,u 4 ) = for u x + 2u 2 + 2u 3 + 3w 4 > s 

(4.12) ^ i k (s,u h u 2 ,u s ,u 4 ) = i k {s,u 1 ,u 2 ). 

«3,li4>0 
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Furthermore we have the recursions 

(u 3 + l)i k (s + 1, ui, u 2 , u 3 + 1, u 4 ) = 

2u 1 Z fc (s - 1,1*!, U 2 , U 3 , M 4 ) 

+ A(u 2 + i)i k (s - 1, iti, u 2 + 1, w 3 , u 4 ) 
+ 4(m 2 + - 1, Mi, M 2 + 1, u 3 - 1, M 4 ) 

+ 4(m 2 + - 2, Mi, M 2 + 1, U 3 - 1, M 4 ) 

+ 2(u 3 + l)z fc (s, Wl, "2, ^3 + 1, U 4 ) 
+ 2u 3 i k (s - 1, Ml, M 2 , M 3 , M 4 ) 

+ 6(u 3 + l)i k (s - 1, m, u 2 , u 3 + l, u 4 ) 
+ 2(m 3 + l)i k (s - 2, Mi, m 2 , m 3 + 1, M 4 ) 
+ 2n 3 i fc (s - 2, ui, u 2 , u 3 , u 4 ) 
+ 4(u 4 + l)ik(s, Wi, u 2 , u 3 - 1, u 4 + 1) 
+ 4(m 4 + l)i k (s - 1, Mi, m 2 , m 3 - 1, m 4 + 1) 
+ Au 4 i k (s - 1, ui, u 2 , u 3 , u 4 ) 

+ 4(M 4 + l)i k (s - 1, Mi, M 2 , M 3 , M 4 + 1) 

+ 4u 4 i k (s - 2, ui, u 2 , it 3 , u A ) 

+ 2(m 4 + l)i k (s - 2, Mi, M 2 , M 3 , M 4 + 1) 

+ (2s - 2mi - 4m 2 - 4m 3 - Qu 4 )i k (s, m x , m 2 , m 3 , m 4 ) 

+ 2(2(s - 1) - 2ui - 4u 2 - 4u 3 - Qu 4 )i k (s - 1, ui, u 2 , u 3 , u 4 ) 

(4.13) + (2(s - 2) - 4m 2 - 4m 3 - 6u 4 )i k (s - 2, m x , m 2 , m 3 , m 4 ) 

and 

2(m 4 + l)i k (s + 1, Mi, m 2 , m 3 , m 4 + 1) = (m 3 + l)ik(s, Ui, m 2 , m 3 + 1, m 4 ) 

(4.14) + 2(m 2 + l) k (s, ui, u 2 + 1, u 3 , u 4 ). 



The proof of Lemma H] is presented the SM and follows by removing a specific arc in a 
labelled C 3 -element and careful accounting of the resulting arc-configurations. 

Proposition [5] and Lemma H] put us in position to compute the generating function of 
colored V fc -shapes 

(4.15) I k {x,y,z,w,t)= ik(s,u)x s y ul z U2 w u H Ui . 

S,1£l, 1*2,1*3, 1X4 
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Proposition 6. Fork > 2, the generating function of colored shapes is given by 

, 1 + x (xil + (2w - l)x + (t - l)x 2 ) 
(4.16) lk( x , V, z, w , t) = — - — F fc 



\ 8 2 
where 9 = 1 — (y — 2)x + (2w — z — l)x 2 + (2w — z — l)x 3 . 

Proof. The first recursion of Lemma H] implies the partial differential equation 

^ = ^i(2x 2 + 4* 3 + 2x 4 ) - ^(2xy + 2* 2 y) 
aw ax ay 

dl 

+ —^{-Axz + Ax 2 w + 4x 2 - 4x 3 2; - 8x 2 z + 4x 3 w) 
dl 

+ — — (— 4xw + 2x — Qx 2 w + 6x 2 — 2x 3 w + 2x 3 ) 
dl 

(4.17) + — -(-6xt + Axw - 8x 2 t + Ax 2 w + 4x 2 - 2x 3 t + 2x 3 ). 

at 

Analogously, the second recursion of Lemma|4]gives rise to the partial differential equation 

dlk dlh dlh 

(4 ' 18 > 2 af = ^ + a7 21 ' 

Aside from being a solution of eq. (I4.17P and eq. (14.181) . we take note of the fact that 
eq. (I4.12p of Lemma H] is equivalent to 

(4.19) l k (x,y,z,l,l) = W k (x,y,z). 

We next show 

• The function 

„/ ,n 1 + x 

y, z,w,t) - x _ ^ _ ^ x +{2w _ z _ ^ + {2w _ z _ 1)a , 3 x 



x(l + (2w - l)x + (t- l)x 



2\ 



\l-(y- 2)x + (2w-z- l)x 2 + (2w-z- l)x 3 ) 2 , 

is a solution of eq. (I4.17P and eq. (I4.18p . 

• its coefficients, il(s, u±, «2, M3, M4) = [x s ?/ ul 2; M2 w n3 t M4 ]I^(x, y, z, w, t), satisfy 

il(s, U\, « 2 , «3, 1*4) = for Ux + 2u 2 + 2u 3 + 3u 4 > s, 

• = W fc (x,y,2;). 

We verify by direct computation that l* k (x, y, z, w, t) satisfies eq. (I4.17P as well as eq. (I4.18p . 
Next we prove i k (s, u±, 112, M3, W4) = for u\ + 2w 2 + 2w 3 + 3«4 > s. Since I%(x, y, z, w, t) is 
analytic in (0, 0, 0, 0, 0), it is a power series. As the indeterminants y, z, w and t appear 
only in form of products xy, x 2 z or x 3 z, x 2 w or x 3 w, and x 3 t, respectively, the assertion 
follows. 
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Claim. 

IfcO, V, z, w, t) = I fc (x, y, z, w, t). 
By construction, i* k (s,u) satisfies the recursions eq. (I4.13P and eq. (14. 14ft as well as 
il(s, Ux,U2, u 3 , 114) = for u\ + 2u 2 + 2«3 + 3«4 > s. Eq. (14. 19}) implies 

/] i* k {s, Ui, u 2 , u 3 , w 4 ) = ik{s, Mi, m 2 ). 

«3,1X4>0 

Using these propertites we can show via induction over s 

V s,u 1 ,u 2 ,u 3 ,u 4 > 0; it(s,u 1 ,u 2 ,u 3 ,U4) = i k (s,u 1 ,u 2 ,u 3 ,u A ) 
and the proposition is proved. □ 

5. The main theorem 

We are now in position to compute Q k (z). All technicalities aside, we already intro- 
duced the main the strategy in the proof of Proposition HI as in the case k = 2 we shall 
take care of all "critical" arcs by specific inflations. 

Theorem 3. Suppose k > 2, then 

(5.1) Q*(*) = fL±^ F fc (0(g)) , 

where 

q{z) = I - z - z 2 + z 3 + 2z* + z 6 - z 8 + z 10 - z 12 

z\l-z 2 -z' + 2z"-z 8 ) 

5.2 ${z) = — . 

q{z) 2 

Furthermore, for 3 < k < 9, Qfc(n) satisfies 

(5.3) Q fc (n) ~ c fc ra -((^-i) 2 +(fc-i)/2) 7 -™ ? ; or some Cfc > 0; 

where 7^ zs t/je minimal, positive real solution of i9(z) = p\, see Tabled 



k 


3 


4 


5 


6 


7 


8 


9 


6{n) 


n~ 5 


'21 

n 2 


n -18 


bb 

n 2 


n -39 


lUb 

n 2 






2.5410 


3.0132 


3.3974 


3.7319 


4.0327 


4.3087 


4.5654 



TABLE 2. Exponential growth rates j k and subexponential factors 6(n), for 
modular, /c-noncrossing diagrams with minimum arc length four. 



Proof. Let denote the set of modular, /c-noncrossing diagrams and let 3 k and Jk(s,u) 
denote the set of all V^-shapes and those having s arcs and Ui elements belonging to class 
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Cj, where 1 < i < 4. Then we have the surjective map, 

inducing the partition = U 7 </?fc (7), where ip~j~ (7) is the preimage set of shape 7 under 
the map (fk- This partition allows us to organize Qfc(#) with respect to colored V^-shapes, 
7, as follows: 

(5.4) Qk(z) = J2 E ^(4 

We proceed by computing the generating function Q 7 (z) following the strategy of Propo- 
sition HI also using the notation therein. The key point is that the inflation-procedures 
are specific to the Cj-classes. In the following we will inflate all "critical" arcs, i.e. arcs 
that require the insertion of additional isolated vertices in order to satisfy the minimum 
arc length condition. In the following we refer to a stem different from a 2-stack as a 
t-stem. Accordingly, the combinatorial class of t-stems is given by (M — Jl 2 ). 

• Ci-class: here we insert isolated vertices, see Fig. [7J and obtain immediately 

Fig. 7. Ci-class: insertion of at least three vertices (red) 

(5.5) d(z) = 

• C 2 -class: any such element is a pair ((i, i+2), i + 3)) and we shall distinguish 
the following scenarios: 

— both arcs are inflated to stacks of length two, see Fig. [HI Ruling out the cases 
where no isolated vertex is inserted and the two scenarios, where there is no 
insertion into the interval [z + 1, i+2] and only in either [i, i+ 1] or [i + 2, i + 3] , 
see Fig. [HJ we arrive at 

C^ a) = IR 4 x [(Seq(Z)) 3 - £ - 2(Z x Seq(Z))]. 

This combinatorial class has the generating function 

^-'(fe)'- 1 -!^)- 

— one arc, (i + 1, i + 3) or (i, i + 2) is inflated to a 2-stack, while its counterpart 
is inflated to an arbitrary f-stem, see Fig. [9] Ruling out the cases where no 
vertex is inserted in [i + 1, i + 2] and [i + 2, i + 3] or [i, i + and [i + 2, i + 3], 
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ABC 



Fig. 8. C2-class: inflation of both arcs to 2-stacks. Inflated arcs are colored 
red while the original arcs of the shape are colored black. We set A = [i + 1, i + 2] , 
B = [i+2, i+3] and C = [i+2, i+3] and illustrate the "bad" insertion scenarios as 
follows: an insertion of some isolated vertices is represented by a yellow segment 
and no insertion by a black segment. See the text for details. 

we obtain 

C 2 b) = 2ft 2 x (M - ft 2 ) x ((Seq(Z)) 2 - £) x Seq(Z), 



having the generating function 




Fig. 9. C2-class: inflation of only one arc to a 2-stack. Arc-coloring and labels 
as in Fig. [8] 

— both arcs are inflated to an arbitrary f-stem, respectively, see Fig. [lOj In this 
case the insertion of isolated vertices is arbitrary, whence 

C<, c) = (M-ft 2 ) 2 x (Seq(Z)) 3 , 

with generating function 




ABC 



Fig. 10. C2-class: inflation of both arcs to an arbitrary f-stem. Arc-coloring 
and labels as in Fig. [8] 



As the above scenarios are mutually exclusive, the generating function of the C2- 
class is given by 

C 2 {z) = C^\z) + cf\z) + Cf{z). 

Furthermore note that both arcs of the C2-class are inflated in the cases (a), (b) 
and (c). 

• C^-class: this class consists of arc-pairs (a, (3) where a is the unique 2-arc crossing 
and has length at least three. Without loss of generality we can restrict our 
analysis to the case ((i, i + 2), (i + 1, j)), (J > i + 3). 

— the arc (i + l,j) is inflated to a 2-stack. Then we have to insert at least one 
isolated vertex in either + 1] or [i + 1, % + 2], see Fig. [TTJ Therefore we 
have 

d, a) = R 2 x (Seq(Z) 2 - £), 
with generating function 

Note that the arc (i, i + 2) is not considered here, it can be inflated without 
any restrictions. 

— the arc (i + 1, j) is inflated to an arbitrary t-stem, see Fig. [TT]) . Then 

e^ b) = (M - R 2 ) x Seq(Z) 2 , 
with generating function 




A B 

Fig. 11. C3-class: only one arc is inflated here and its inflation distinguishes 
two subcases. Arc-coloring as in Fig. [8] 

Consequently, this inflation process leads to a generating function 

C 3 (z) = Ct\z) + C^(z). 
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Note that during inflation (a) and (b) only one of the two arcs of an C3-class 
element is being inflated. 
• C 4 -class: this class consists of arc-triples (ai,/3,a 2 ), where a,\ and a 2 are 2-arcs, 
respectively, that cross 0. 

— (3 is inflated to a 2-stack, see Fig. [12j Using similar arguments as in the case 
of C3-class, we arrive at 

ej a) = ft 2 x (Seq(Z) 2 - £) x (Seq(Z) 2 - £), 

with generating function 

— the arc is inflated to an arbitary f-stem, see Fig. [121 

e< b) = (M - ft 2 ) x Seq(X) 4 , 
with generating function 




ABC D 

Fig. 12. C4-class: as for the inflation of C3 only the non 2-arc is inflated, 
distinguishing two subcases. Arc-coloring as in Fig. [8] 



Accordingly we arrive at 
(5.8) C 4 (^) = Ci a) (^ + Ci b) (^). 

The inflation of any arc of 7 not considered in the previous steps follows the logic of 
Proposition HJ We observe that (s — 2u 2 — u 3 — u 4 ) arcs of the shape 7 have not been 
considered. Furthermore, (2s + 1 — u\ — 3^2 — 2^3 — 4u 4 ) intervals were not considered for 
the insertion of isolated vertices. The inflation of these along the lines of Proposition H] 
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gives rise to the class 

§ = J^s-2u 2 -u 3 - Ui x ( SEQ (£))2.+l-ui-3u a -2u3-4u4 j 

having the generating function 

s—2u2 —U3 —U4 



S(z) = 1 — x 2 — ^ 1 

1 ~ T=? (fe + (1=5) 



\ 2s+l — u\— 3it2— 2113— 4«4 



1 - 

Since all these inflations can freely be combined, we have 

q 7 = e? 1 x q u 2 2 x e^ 3 x x s, 



whence 



Q 7 (*) = Cx^r • C2W 2 • C 3 (z) U3 • C 4 W U4 • S(*) 

1 



1 — z 



where 



z 4 



'° {Z) - l-2z + 2zz-z±-2zs + z^ qi{z) ~ Z 

, , 2(1 - 4z 3 + 2z 4 + 8^ 5 - 6z 6 - 7z 7 + 8^ 8 + 2z 9 - 4z 10 + z 11 ) 

ft (2) = : 

1 — Z 

C 3 (*) = z{2 - 2z 2 + z 3 + 2z 4 - z 5 ) 
q(» = z 2 (5 - 4z - 3z 2 + 6z 3 + 2z 4 - 4z 5 + z 6 ). 
Observing that Q 71 (z) = Q 72 (z) for any 71,72 G Jk(s,u), we have according to eq. (15. 4ft : 

s,u>0 

where u > denotes it, > for 1 < i < 4. Proposition [6] guarantees 

ik(s,u) x n y Ul z U2 w u H u * 

s,u>0 

1 + X 



x 



1 - (y - 2)x + (2w - z - l)x 2 + (2w-z - l)x 3 
x(l + (2w - l)x + (t - l)x 2 ) 
(1 - (y - 2)x +(2w-z- l)x 2 + (2w-z- l)x 3 ) 2 
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Setting x = <sq(z), y = qi(z), z = ^(z), w = q 3 (z),t = q^(z), we arrive at 

Q , , l-z 2 + z* 

^ k[Z) l-z-z 2 + z 3 + 2z* + z G -z s + z w -z 12 X 

z\l -z 2 -z i + 2z 6 - z s ) 



{l-z-z 2 + z 3 + 2z i + z e -z 8 + z 10 - z 12 ) 2 , 

By Lemma [H Qfc(-s) is D-finite. Pringsheim's Theorem [17] guarantees that Qfc(z) has 
a dominant real positive singularity 7^. We verify that for 3 < k < 9, 7^ which is the 
unique solution of minimum modulus of the equation d(z) = p\ is the unique dominant 
singularity of Qk(z), and ~&'(z) 7^ 0, see the SM. According to Proposition [2] we therefore 
have 

Q fc (n) ~ c k n -((fc-i) 2 +(fc-i)/2) ( 7 -i)n j for some Cfc > 
and the proof of Theorem [3] is complete. □ 

Remark 1. We remark that Theorem [3] does not hold for k = 2, i.e. we cannot compute 
the generating function (^2(2) via eq. (15. ip . The reason is that Lemma H] only holds for 
k > 2 and indeed we find 

1-^ + ^ /z 4 (l - z 2 - z* + 2z 6 - z 8 ) 
(5.9) Q 2 (z) ^ r~\ — F 2 



g(z) V l( z ) 2 

However, the computation of the generating function (^2(2) in Proposition H] is based on 
Lemma [21 which does hold for k = 2. 
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